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LINEAR DIFFERENCE EQUATIONS AND THEIR ANALYTIC 

SOLUTIONS* 

BY 

R. D. CARMICHAEL 

Introduction. 

The first existence theorem for difference equations in which analytic solutions 
were treated was obtained by Guichard f who, in a paper published in 1887, 
proved that if <f>(x) is any entire function whatever there exists another entire 
function f(x) satisfying the relation 

f(x + l)-f(x) = <t>(x). 

Demonstrations of this theorem have also been given by Appell J and Huewitz. § 
Hurwitz further showed how to solve the equation 

4>(x)f(x + l)- X {x)f{x) = f{x), 

where <f>(x), x( x )< ^(x) are meromorphic functions. 
Mellin II has considered the nature of the solution of 

f(x + l)-r(x)f(x) = s(x), 

where r ( x ) and s ( x ) are rational functions. 

Barnes ^[ has contributed an interesting study of the equation 

<f>(x)f(x+<o)- X (x)f(x) = yfr (x), 

where m is any constant different from zero and 4>(x), %(#), yfr(x) are entire 
functions. 

Watson** has considered the equation of the second order 

A(x)f(x+ 1 ) - B(x)f(x) + C(x)f(x - 1) = 0, 



* Presented to the Society September 6, 1910. 

tAnnales de l'Eoole Normale superieure, ser. 3, vol. 4 (1887), p. 361. 
% Journal de Matheroatiques, ser. 4, vol. 7 (1891), p. 157, especially chapter I, pp. 
159-176. 

§ Acta Mathematica, vol. 20 (1897), p. 285. 

|| Acta Mathematica, vol. 15 (1891), p. 317. 

If Proceedings of the London Mathematical Society, ser. 2, vol. 2 (1905), p. 438. 

**Proceedings of the London Mathematical Society, ser. 2, vol. 8 (1910), p. 125. 
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with some restrictions on A(x), B(x), C(x), the object of his investigation 
being to obtain a solution which is of use in applications to the theory of linear 
differential equations. Two analytic solutions are obtained, but the question of 
their independence is left unanswered. 

Finally, Galbeun * has employed the Laplace transformation which carries a 
linear difference equation with polynomial coefficients into a differential equa- 
tion, and by means of this transformation has obtained important existence 
theorems for the case in which the difference equation has polynomial coefficients. 

The primary object of the present paper is to prove that there exists a funda- 
mental system of simple analytic solutions of the system of n linear homogeneous 
difference equations 

n 

(A) G i (x + l) = x a J^a i .(x)G j (x) (t-=l, ••-,»), 

.7=1 

the known quantities entering into the equations being defined as follows : 

1) a is a constant ; 

2) the functions a r (x) are single- valued and 

a iJ (x) = a iJ + a' ij x- 1 + a'/.x- 2 •■■ (i, j = 1, • • •, n) |x|>is?'; 

3) the constants a., are such that the roots of the equation 

1 v « n ' « I 1 if % =J, 

are in absolute value different from each other and from zero. 

In § 1 a form of the method of successive approximation is developed in con- 
nection with a single equation of the first order and by means of it a study is 
made of the nature of a simple particular solution of the equation. 

The methods of § 1 are readily extended in § 2 in such a way that they lead 
to a similar solution of the normal system 



(B) f i (x + l)-a i (l+-) f i (x) = T,U*)/A*) 

< I a l | < | a 2 1 < • • • < | a n | 



(*=!.-■ 



where 

and 

<l> ij .(x) = <t>' i ' J x- 2 +<f>' i j'x- 3 + ... (l,j = l, .-.,n) |aj|^i?>l. 

A method of induction in § 3 leads to the determination of a fundamental 
system of solutions of (B) and to the development of important properties of 
the functions of which it is composed. 

In § 4 it is shown that system ( A ) can be reduced to the normal form ( B ) 

*Comptes rendus, vol. 148 (1909), p. 905; vol. 149 (1909), p. 1046; vol. 150 (1910), 
p. 206; vol. 151 (1910), p. 1114. 
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and thus the general existence theorem for (A) is obtained. Finally, an appli- 
cation of this result is made in deriving the general existence theorem for a 
single equation of the nth order. 

I desire to express my cordial thanks to Professor G. D. Birkhoff for valuable 
suggestions and encouragement in the preparation of the paper. 

§ 1. The method of successive approximation for a single homogeneous 
equation of the first order. 

We shall begin by considering a single homogeneous equation of the first 
order in the normal form * 



(1) 



or 



f(x + l)-a(l+ ljf(x) = *(*)/(«), a * 0, 



f(x + l) = ir{x)f(x), +(x) = a(l+iy+4,(x), 

where 

(2) <£(cc) = tf>" x- 2 + <t>'"x-' i + ..., for \x\^B>l. 

We assume that ^r(x) is single- valued and analytic at all points of the complex 
plane where it is defined, and that that determination of (l + l/a;) m is taken 
which approaches 1 asa; approaches infinity. The sequence of equations 

fV(x + l)-a(l + iyfn(x)~0, 

/w (a . + 1) _ a (l + 1 )">>(*) = 4>W m (*), 

f w (x + 1) _ a (l + ljf^(x) = <Kx)J™{*)> 



(3) 



will serve to define a sequence of approximation f unctions / ( ' ) ( a; )>./' <2) ( x )i 

A particular solution of the first equation is evidently a x x m . All the other 

equations of the set have the general non-homogeneous form 

(4) g( x+ \)-a(\+~) g{x) = v{x). 

* It is easy to reduce the general first order equation 

G ( x + 1 ) = z» ( B + f/z-i +6"x-i H ) G {x ) , e + o, 

to this form. Put 

G(x) = {T(x)} a f(x). 

Then 

Hz + i) = (o + e'*-i + o"x-* + - ••)/(*); 

and for \x\ sufficiently large this evidently may be put in the form (1). 
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One readily obtains a formal solution of this equation as follows : Let 

ff(x) = a x x m g(x), 

where a x x m has any determination fixed by making a cut along the negative axis 
of reals in the 35-plane. Substituting in (4) this value of g{x) and dividing by 
a x+1 (x + l) m one has 

g(x + l)-g(x) = a , +l V { ( x X +1)m , 

an equation which has the formal solution 

_, .__ V (x ) v( x + 1) v(® + 2) 

9K X )- a x+1 (x+l) m a x+2 (x+2)»' a x+3 (x + 3) m "" 

Hence 

77(35) 77(35 + 1) 77(33 + 2) 

(5) 9^= 



Now 



where 



■Nr -i^iy «-K)" 

r M+-j =\(x)\(x + 1) ... \(x +r — 1), 



and therefore (5) may be written in the form * 

, fi v ( x_ v ^±i) 

W fl'W- ^ \(x)\(3j + l)... \{x+i)' 

This yields not merely a formal but an actual solution for any values of 35 for 
which it converges ; and if 77 ( x ) is analytic, then (by Weierstrass's theorem) it 
yields an analytic solution throughout any two-dimensional continuum where it 
converges uniformly. 

By the aid of (6) we obtain a sequence of formal solutions f {t> {x) satisfying 
(3) identically. Employing the notation 

S(x)=-± ^^ 

and using for each equation that particular formal solution which is obtained by 
adding the solution a x x m of the homogeneous equation to the solution derived 
from (6), we have 

f<- l \x) = a x x m , 

f<n ( x ) = a x x m + S x ( 4> a x x m ) , 

fW(x) = a x x m + SJ<f>a x x m + <f>S x (4>a x x m )} , 



*Cf. Mellin, Acta Mathematica, vol. 15 (1891), p. 319. 
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It is easy to see that this leads to a formal solution f x ( x ) of (1), denned by the 
series 



(7) 



f l (x) = a*x m + S x (<f>a*x m )+ S x {<t>S x (<f>a*x m )} 



+ SA4>S x {cj>S x (<f>a^)}] + .... 

That this series formally satisfies (1) we may verify by forming the difference 
/*j(a;+ 1) — a(l + l/x)" ! i /* 1 (a;) and effecting a reduction of the resulting series 
by means of the general formula 



s x+l (x)-a{i + ljs x ( x ) = x- 



Consequently whenever all the series S x in the terms of (7) and also the series 
(7) itself converge, this series represents a solution of (1). 

In studying the nature of the formal solution (7) (and in the work of the 
succeeding sections also) it will be convenient to have names assigned to certain 
regions of the plane about which we shall have to speak often. Let R be the 




Fig. 1. 

circle | x | = R. Let I be a line cutting the positive axis of reals at right angles 
and not cutting R. The part of the finite plane to the right of and including 
I will be called a jP-region. Let m and n be parallel to the axis of reals, m 
lying above and n below R. The part of the finite plane lying above and 
including m will be named a ^-region ; that lying below and including n, a 
Q-region ; these two together, a (^-region. A region formed by the combina- 
tion of the R- and (^-regions will be called a jD-region. The part of the finite 
plane not included within D will be called the -E'-region corresponding to D. 
It is scarcely necessary to remark that Q and Q have each a part in common 
with P. 

In what follows immediately we are to prove the existence ofy^a;) and develop 
its properties in particular regions R l , Q l , D x , i?, ; for convenience all these 
regions are defined at once. For any positive constant c , less than | a | , there 
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exists an R' such that 

c , for I x I > S ' . 
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From the nature of the function <\>(x) it follows that there exists an M such 
that 

Mc 
[ d> ( x ) I < i — r , , f or I x I S i? . 

\t\ j\ j as | 2 * ii 

If one writes 

a; = m ± ry/— 1 

where « is positive and both u and -y are real, it is now clear that one may choose 
particular regions P, Q, D, E — call them P l , Q lf D l , E l respectively — such 
that 

Mc / „ 1 V 



(8) 

(9) 
(10) 



\4>{x)\ 



x 



c , for x iu D l ; 



A 

v 



Mm 



+ 



2Jf 



>Jf+l, for a: in P t ; 

2 < 1 , for x in Q t , A being a constant ; 



(11) a(l+j +<f>(x)j=0, for x in Z^. 

Let us consider the function 

(£(« + i)a* +i (a; + i) m 



(12) £ (<£a*af) = y; , , , w 

V ; xVr ; ^J X(sB)X(as+l)...\(aj+») 

Employing inequalities (8) one has the relation 



= £ 



<j>(x + i)a x x m 
\(x + i) 



(13) 



<£(a; + i)a*x m 



\(x + i) 



«£ 



M\a x x m \ 



subsisting throughout D x . The series of the last member of (13) is evidently 
uniformly convergent when x lies in any closed region in D x . The same then is 
true of the series in the first member, and hence also of the series in (12). 
Therefore S x (<f>a x x m ) is a uniformly convergent series in the same closed region 
and hence, in accordance with a well-known theorem of Weierstrass, represents 
in this region an analytic function of x . From this it follows that the first and 
second terms of the series in (7) represent functions which are analytic in D x . 

In order to show that the remaining terms are analytic in the same region it 
is necessary to know a bound to the magnitude of \S x (<f>a x x m ) | as x varies. In 
carrying out the needful reckoning it is convenient to treat the regions P x and 
Q x separately. 

* The symbol <^ is to be read is less term by term than. 
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First, let us consider the region P x . Of the two series 

00 M °° M 
y. -, TT? -, , T2 t ^r? 5 f° r * in P. , 

the terms of the first in order are equal to or less than the corresponding terms 
of the second. But the latter series is less than the integral 



r 



M M 

, M _!Z 2 «-l 

and therefore 

M 

(14) S x (<f>a x x m )<\a x x m \^-^. 

Hence 

M 

I S(<j)a x x m )\<\a x x m \ ^ . 

Using this result and making an argument similar to that associated with 
inequality (13) one may show that S x {<f>S x (<f>a x x m )} is analytic throughout 
P x . Again, in view of (14) we have 

- - - f i| 1 . M 2 

Similarly it may be proved that S x \_4>S x {$S x {4>a x x m )}'\ is analytic in P l ; 
an upper bound for the value of S x [<£# x ( <f>/S x (<f>a x x m )}~\ may be obtained; 
and so on. Hence in series (7) every term is analytic throughout P 1 ; more- 
over the following inequality has been established for x in JP t : 

S x (<}>a x x m ) + S x {<t>S x (<f>a*x™)} + S x [<}>SJ<l>S x (<t>a x x m )}-\+ ■■- 

( 15 ) ^M\a x x m \ M 2 \a x x m \ M s \a x x m \ 

x u-i + (u-if + Ju^rf + " " 

From (9) it follows that the second member of (15) is absolutely and uni- 
formly convergent throughout any closed region in .P, ; and therefore the same 
is true of the first member. Hence series (7) is uniformly and absolutely con- 
vergent throughout this region. But its terms are analytic functions of x ; and 
therefore, by Weierstrass's theorem, its sum /^(x) is an analytic function of x 
throughout P l . Evidently, then, f(x) =f 1 (x) is a solution of (1). 

A comparison of (7) and (15) shows that the following inequalities are valid 
for x in P 1 : 

[MM 2 1 

l^ M l| 1 -,^i-(^i7----!<I/ 1 ^)l 

< a x x m 1 + t + , ^ + • • 

1 [ u — 1 ( u — 1 ) 



,^ 2 

<\a*x m \— 2 . 
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Hence we have 

( 16 ) \™\ VJZlt <\M*)\<\**\^=i> 

a set of fundamental inequalities bounding the increase or decrease of f^(x) 
throughout _Pj . 

Again using relation (13) as a point of departure, we obtain similar results 
valid in Q l . In this region v is different from zero. It is evident that 

f> M » M <0 f p dz M\_Mtt 2M A 

h {u+iy+v^^ ( M+ ^+« 2= 1 M X z^-w + * J ~ "T" + "?" <¥ ' 

J. being a properly chosen constant. Comparing this result with (12) and (13) 
one has 

(17) | S x (<j>a x x m ) | ^ S x (4><far) < | a x af* 1 — . 

On account of (10) and (18) it is clear that one may employ an argument 
similar to that associated with (13) to show that S x {<f>S x (<j>a x x m )} is an ana- 
lytic function of a? in Q v . Again, on account of (17) one has 

(18) SAWM**)} <s\4>\*'* t \*}=* s^tt) 

In like manner it may be shown that the function S x [_<j>S x { <f>S x (<j>a x x m ) } ] 
is analytic throughout Q x ; also, an upper bound for the value of 

S x {<$>S x {<l>S x ($a?x™)}-] 

may be obtained ; and so on. Hence every term in series (7) is analytic in Q 1 . 
Moreover, for this region the following relation exists : 

(19) S x {<$>a x x™)+S x {4>S x {<t>a x x™)} + S x [<j>S x {<j>S x (<l> a *x<«)}] + ••• 

« \a x x m - + a x x m -j + \a x x m -r- + • • • . 

By a method similar to that employed in the previous case it is now easy to 
show thaty^aj) is analytic in Q l , and that throughout this region it is bounded 
as in the inequalities 

( 20) |^ ar |L=^ < | /i(aj) | < | rf . aj »|_^_. 

It remains to determine the nature oif^x) in the region E . For this study 
series (7) does not furnish a convenient means. We proceed as follows : If we 
solve (1) for_/(x) and in the result write x — 1 for x we have 

(21) f(x-l) = f(x-l)f(x) 
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where 

(22) f(x) = {f(x)}-\ 

By means of (21), the solution f x ( x) which is known to exist in P l can be 
extended across the plane to the left ; and in this way the nature ofy^a;) in E y 
is determined. 

Consider the set of points [$] which is formed of the set [a] and of the 
sets [a — i ] , i = 1 , 2 , 3 , •••, where [ a ] is the set of singularities of ty ( x ) . 
A point a — i is said to be congruent to a on the left ; in a similar way con- 
gruence on the right is defined. It is easy to see that /"j (a;) is analytic at every 
point x not belonging to [#]. For, by (21), 

f x {x) = ty(x')ty(x + 1) ... $(x + v— l)f l (x' + v). 

v can be so chosen that x + v lies in P x . For such a value of p, f l (x)is 
analytic at x + p; also, ^jr(x) is analytic at each of the points x, x'+ 1, • • •, 
x + v — 1 . Therefore f x ( x ) is analytic at x. Combining this with previous 
results, we see that f x ( x ) is analytic throughout the finite plane except at the 
singularities [a] of yfr(x) and the points- congruent to them on the left; at 
these points its character depends upon the nature of the particular yjr(x) 
involved. 

Let each point of [&] be inclosed by a small circle of radius e having the 
point as its center. Denote by E u the whole of that part of E x no point of 
which is within one of these small circles. 

Inequalities bounding the increase or decrease of J\ ( x ) for x in E u may be 

found as follows: Consider any positive constants a_ and a + satisfying the 

relations 

«_ < | a | < a + • 
Since 

i{ a ( 1+ ^) + < K a, )}= a ' 

it is clear that r may be chosen so that 



«U> 



:(l + ^ + ^(, ) 



> a_ , for Ice I > r , 



and also so that the circle | x | = r cuts the boundary of E t twice on the left of 
the origin and at no other point. Since ^(x) is everywhere analytic in E l€ it 
follows that there exist constants JVj and N 2 such that 

iVj = \fi(x" ) | < W 2 , for x" in E u and r= \x"\ =r + 1 . 
But (21) yields the equation 

f^x" — fj.) = yfr(x" — (i)yjr(x" — fi + 1) • • • i/r(x" — l)/j(a)"). 
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Hence we have 

iV, a~» S \f(x" — fi) | < W 2 aZ» , for x" in E u and r S | x" \ = r + 1 . 

As fi varies over the range /* = , 1 , 2 , ••-, x"— /j runs over all points in i? 1€ 

and outside of the circle \x\ =r. Also, yj ( x ) is finite in the part of E u 

which is in the circle | x | = r . Hence it is clear that there exist constants iV, 

and N 2 such that 

(23) N l \a* + \S\fi(x)\<N 2 \a*_\, for x in ^ , 

the result sought. 

So far we have said nothing of the position of the zeros of/j(a;). By means 
of equation (1) we shall show that they can occur only at the singularities of 
yfr ( x ) and at points congruent to them on the left. We have 

f x {x' + /*) = f(x' + fi -l) + (x' + M - 2) • • • ^(x')fi(x'), 

where p is a positive integer. From inequality (16) it follows that for every x 
there exists a fi such that f x {x + fi) is different from zero. Hence /,(aj') can 
be zero only when one at least of the factors i/r ( as' ) , ijr(x+l), • • • , i|r(x'+/*— 1) 
is infinite : and therefore the zeros of f t (x) must all be contained in the set of 
points composed of the singularities of yjr(x) and the points congruent to them 
on the left. 

The principal results contained in the section may be enunciated as follows : 

The difference equation 

f(x+l) = f(x)f(x) orf(x-l) = jr(x-l)f(x), +(x) = {ylr(x)}-\ 
in which ^jr(x) is single-valued and 

f(x) = a(l + -\ +<f>(x), <f>(x)=<f>"x- 2 +<t>'"x- 3 +---, for\x\^R>l, 

has a particular solution f(x) which is analytic throughout the finite plane 
except at the singularities [a] of-^r^x) and the points congruent to them on 
the left, and which can vanish only at the singularities [/3] ofyfr(x) and at 
points congruent to them on the left. If D l is defined by relations (8) to (11), 
then throughout any closed region in E x this solution is representable as a 
uniformly converging series in the form 

fi(x) = a*x m + S x (<j>a*x m ) + S x {<f>S x (cl>a*x m )} + S x [<j>S x {<t>S x (<f>a*x m y^ + • . .. 

The increase or decrease of f(x) for x in the corresponding regions P x and 
Q x is bounded as in the inequalities 

u - 2M- 1 , „ , x , , w - 1 

l^l v^rMzri < \A( X )\ <\ a * xm \u _ M -i ' J° r xmP ^ 

v — 2A v 

K^'I^X'^ |/i0»)l <\<?**\^ZrA> for x in Q x . 
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Also, for any positive constants a_ and a + satisfying the relations 

a_ < | a | < a + 
there exist other constants iVj, N 2 such that 

^"iK+l = |/i(*)| <^.K|. for * "» E u . 

Remark. About each statement in the conclusion of this theorem it should 
be observed that, insofar as the statement refers to the region D x or any part of 
it, it is independent of the nature of <f>(x) in the region exterior to D l . More 
explicitly: If<f>(x) is any function which is analytic in some D-region and if 
also 

Mc 
100)1 <uT2' far x in I), 

then there exist regions P x , Q x , D x having as before the properties (8) to (10) 
and there is a particular solution f(x) of the difference equation having in D x 
the same properties as the solution f(x) of the equation in the theorem. These 
conclusions follow from the preceding discussion without change of argument. 

§ 2. The method of successive approximation for a normal system of n equa- 
tions. Properties of the solution thus obtained. 

Let us consider the system of n linear difference equations of the first order 
in the normal form 

or 

JL / l\ m ' 

/ 1 (* + l)-Z^(«)/,(*). *«,(*) aV^ + a) + <M a! )' 
where 

(25) 0<|a 1 |<|a 2 |<-.. <|a.| 
and 

(26) <j> ij (x) = <t>" J x- 2 +<f>':' J x- 3 +... (i,j=l, ...,n), for \x\^R>l. 

We assume that ^^{x) is single-valued. If we solve (24) iov f(x) in terms of 
f(x + 1), i = 1, • • •, n, and in the result write x — 1 for x, we obtain 

n 

(27) f i (x-l) = ^^.(x-l)f J (x) (1 = 1, ■■-,»), 

_ J=1 

^(x) having the form 

(28) M*) = °^ 
where A is the determinant 



A^a^l + lJ'+^x) 
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and «..(cc) is the cofactor of the element in thejth column and the ith row of A . 

The object of this section is to obtain a particular solution of (24) and to study 
its properties throughout the plane. 

On account of (25) and (26) it is clear that there exist particular regions P ir 
Q 1 ,D 1 , E l and constants c , M and A such that the following inequalities are true : 

(30) \<f> (sb)|<- t -^, \ a .ll+-\ > c (i,j = l, •-.,«), for as in Z>, ; 



(31) 



/ l\ Ml | / l\ m 'l 

a \ 1+ x) <a \ 1 + «J ! (» = 2,...,n), forasinZ^; 



(32) u > M + 1 , for x in P, , where x = u ± vV — 1 ; 

i Jfir 2Jf , , . „ 

(33) _ < __ + _ < i ) for « in ft ; 

(34) A+0, for|x|SPSP. 

In the study of our system of equations these regions are treated separately. 
It will be convenient to employ the notation 

y( X -f v) ( IN" 1 ' 

^(x)=--S x<(a>)X|(a>+1) ... x<(g + y) . \(-)-«^l+ 5 J ' 

also we write 

w ; = C { a^x mi , 

where C. is any periodic function of x of period 1 and analytic throughout 

A- 

The form of the method of successive approximation employed in the pre- 
ceding section is readily generalized so as to lead to a particular solution of (24). 
Consider the set of systems each of n linear equations, all except the first being 
non-homogeneous : 

f(p(x + l)-a i (l+iy i fP(x) = Q, 

X J j=l (» = l, ■■■,»). 

/ 1 \ m i n 

/c« (a . + i)_ a ^l + _j ff(x) = ^<f> ij (x)ff(x), 

Proceeding as in the foregoing section, one writes ont in order formal solutions 



(t'=l, ••-,«). 
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of these systems, obtaining thus the following set of approximation functions : 

f?{x) = u r 

f?\x) = u i + S, i (±<t> ij u^, 

/fix) = Ui + s„(i: *<,«,) + s*<{ i^s^tt^) }, 

It is easy to verify by actual substitution that a solution of system (24) is 
obtained in the limit functions (if they exist) of the sequences f ( P, ff, f ( ?\ 
• • • , i = 1 , • • • , n; that is, in 

flx)=u i+ 8 xi ( £ 4>.u. ) + si £ ^ ( £ <*>, A ) j 

+ - s -[S^-ll^ A (S te )}] + - 

This solution is in general illusory on account of the divergence of certain series 
contained in the terms of the expansions formally representing the /'s. In 
fact, let a and a be any two a's such that la 1 > I a I . Put 

7 r s *J I r I *^ I s I 

<7 i= =0, i + r; C r =l. 

For i = s the second term in the series in (35) becomes 

„ , N ^,6 (x + v)a* +v (x + v) m <- 

A„W.,«,« >- 2. X (*)... x(x + *) 
(36) 



a* + "(a: + !/)*"• ' 

a series which is evidently divergent at least for some values of <f> sr (x) as for 
instance <f> sr (x) = 1/x 2 . 

On the other hand, there is at least one particular convergent solution, arising 
from the following special choice of values of the periodic functions C t : 

<7 i= 0, j + l; C 1 = l. 

If we call this particular solution f a ( x ) , we have formally from (35) 



f a (x) = S a alx^ + S xi (<t> n a*x™>) + S xt { £ ^(^ a* *•»>)} 
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That all the infinite series contained in (37) are actually convergent in D l is 
readily shown. Consider the series 



(38) 



R ( Mc \ x \Mc" _ / Mc \ 

.11 -"\ nx 2 1 J •" | war fr[ XJ \ nx 2 1 J 



+ 



- V3Ic » - f JkTc ™ - /Jib \n 

" L «* JTi J { nx fr[ "' \ nx' 1 ) J J 

where £ . denotes what S ri becomes when the terms of the series which it indi- 
cates are replaced by their absolute values. On account of the first inequality 
(30) it is evident that all the infinite series in (37) are convergent provided that 
all those in (38) are likewise convergent : for if these expansions are written out 
in detail, the former, term by term, is less in absolute value than the latter. 
Compare (38) with the following series 



(39) 



R / Mc \ - f Mc - ( Mc \\ 

a\ x^ + 8 A \-j a\x^ J + 8 A j -^ 8 A (^ -p a? *"" J J 

- [Mc - (Mc - (Mc \ 1 



+ 



On account of inequality (31) it is clear that the convergence of (39) carries with 
it the convergence of (38) ; this is at once apparent from a term by term com- 
parison of the two series, since (39) is obtained from (38) by replacing each de- 
nominator by one at least as great. But (39) (except for its first term) is of 
the same form as the series in the first members of (15) and (19). Hence, just 
as in the case of the latter series, one proves that throughout the region D x the 
successive terms of (38) represent analytic functions and finally that the series 
itself represents an analytic function. 

We shall now find an upper bound for the values of the n — 1 functions 
f a ( x ) , i = 2 , • • • , n , as x varies within P l ; and, for similar variation of x , both 
upper and lower bounds for the values of f n (x). We have evidently 

- ( Mc \ - f Mc - ( Mc \\ 

fa (») < ^1 {^ a > mi ) + M 1?" Sxl { 1* ^^ ) 1 + " ' 

AM M 2 1 

the latter inequality being deduced by aid of (15). Hence it is clear that 

M 

(40) \fa(K)\<\<*> mi \^TM~\ (» = 2,. .-,«). 

Similarly, one has 

M-l 



|/i.(»)l<K«'*l( 1 + u _^_l ) = l g T ! 



l„_ M-V 
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and also, 

Combining these two inequalities we get 



u-2M-l 



, u ~ 2Jf— 1 , , , ,, , , u — 1 

( 41 ) K^l tc-if-^ l/nHKK^I^rjgin- 

The desired results for the region JP l are contained in inequalities (40) and (41). 
Starting from (19) and working in a similar manner one obtains the follow- 
ing inequalities valid in Q x : 

( 42 ) \fi l (x)\<\at* mi \^A (i = 2, ■■-,»); 

(43) | a\x^ | ^^ < \f n (x)\ < 1 -aJ- 1 ^ ^ :i . 

By means of (27) we are able to determine the nature of f a {x), i = 1 , •••,«, 
in the region E ; for we have 

( 44 ) A(^-i) = i:^(x-i)y;,(x). 

.7=1 

Now the functions /.,( x ) are analytic throughout P x . This solution may there- 
fore be extended across the plane to the left by means of (44) ; it is evident 
from the nature of the process that the singularities of/ a (a;) for x in E x are at 
the singularities of i/r y (x), i,j =1, ■ ■ -,n, and the points congruent to them on 
the left. The nature of the singularities which the functions f a (x) may have 
depends upon the particular functions yjr r (x) entering into a given particular 
system of equations. In special cases it may happen that some or all of the 
functions have no singularities at the points indicated. 

Inequalities bounding the increase or decrease of \f a (x) | for x in E x are not 
as complete as those for x in P x and Q l . We shall however obtain bounds for 
|/ a (a:)| when x is in E l€ ,E u being defined as follows: Let [/S] be the set of 
points composed of all the singularities of i^\.(x), I, j = 1, • • •, n, and the 
points congruent to them on the left. Enclose each point of [/S] by a circle 
of radius e having the point for center. The part of E x no point of which is 
interior to one of these circles is the region E le . 

Eeferring to the definition of ^r {j {x) in equation (28), it is easy to see that 



L a..(x) = 



0ifj=M, 

n 

a^ 1 II a, itj = i, 



Trans. Am. Math. Soc. 8 



114 E. D. CARMICHAEL : DIFFERENCE EQUATIONS [January 

and that 
Hence we have 

From this it follows that for every positive constant a such that a < | a t | there 
exists a positive constant r for which 

n 

2|V r <,( aj )|< a_1 (i = l, ■•-,») |a;|Sr. 

It will be assumed that r is so taken that the circle | x [ = r cuts the boundary 
of E y in two points on the left of the origin and in no other point. 

Since each of the functions f a ( x ) is analytic throughout E u it follows that a 
constant N exists such that 

\f n (x") | <i^for x" in E u and r = \x"\ = r + 1. 

But from (27) we have 

w n n 

f n (x" -?)- z^(*" - /*) E^(*"- /* + 1 )• • • Z ^ v _ UM K-i)4K)- 

Hence we get 

\f il {x"- ( i)\<Na-». 

As /u, varies over the range \i = , 1 , 2 , ■ ■ ■, x" — /u. runs over all points in .Z?,,, 
and outside of the circle |a;| = r. Also, j/^Ca;)! is finite in the part of E u 
which is in the circle |cc| = r. Hence it is clear that there exists a constant 
N, such that 

\f,{x)\<N\a x \ (»=1, ■■■,») asin^.. 

The principal results of the section may be stated as follows : 
The system of difference equations 

n n 

3=1 j=l 

in which yjr..(x) is single-valued and 

+ { .(x) = S i .a i (l + iy i + 4>..(x), 4> i .{x) = & J x-*+W J x->+... 

for\x\=lB>l, 

S.j = when i =(= j ; 8.. = 1 , 

A.as a particular solution, f a (x) , each function of which is analytic throughout 
the finite plane except at the singularities [a] of yjr { .(x) and the points con- 
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gruent to them, on the left. If D l is defined by relations (30) to (34), then 
throughout any closed region in D^ this solution is representable as a uni- 
formly converging series in the form 

f n (x) = S n a^ + S H {4> n a\x^) + 8 xi { ± <£.. S^htf**)} 

1 j= i J 

+A*[i:* tf ^{i:^A(*»«T^)}] + --- «=i, ■•■-*)■ 

:7%e increase or decrease of these functions for x in the corresponding regions 
P l and Q l is bounded as in the inequalities 

u - 2M- 1 , , , s , , . u — 1 , . „ 

K«""'| M _jf_i <l/n( a! )l<l«r a; "' 1 l ! 7Zrjf^l' forxinP^, 

|/ <1 (g)|<jqtai" 1 l tt _ jf _ 1 (i=2,...,n), for x in P^, 

\fn( x )\<\ a i x '" 1 \—-2 (i = 2,...,»), forxinQ^ 

Moreover for any positive constant a such that a < | a t | <Aere exists another 
constant JV such that 

\f. l (x)\<A T \a :e \ (t=l,...,ra), forxinB u . 

Remark. It is important for later use to point out the fact that, in so far as 
any statement in the conclusion of the theorem refers to the region D^ or any 
part of it, it is independent of the nature of <j>. . ( a;) , i,j= 1 , • • • , n , in the region 
exterior to Z) L . Hence, if the functions <j> r (x) are analytic in some D-region 
and satisfy the inequalities 

Mc 

( 45 ) Wn^^nTxl* (*,j = l, ■•-,«), forxinD, 

then there exists a region D x having -properties as in the theorem, and the sys- 
tem of equations has a particular solution f n (x) having in D x the properties of 
f n (x) of the theorem. 

§ 3. Existence of n independent solutions of the normal system. 

For the system of equations considered in the preceding section we shall now 
obtain n particular solutions which are independent and therefore constitute a 
fundamental system of solutions. 
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Since by (37) it is seen that f n (x) is not identically zero, we may make the 
substitution 

(46) gi (x)=Mx), 

(47) ^(*) =/,(«) -^/.(*) (. = 2, ■">«), 

of which the inverse is 

f (x) 

(48) M«) = g l (x), /<{«)= j-j^frW + ffiW (i = 2,--,»), 

where / H (a:), i = 1, •••,«, is the solution of (24) obtained in the preceding 
section. Then we have 

?.(» + !)-/,(»+ 1 )- /° ( J+lj /.fa + 1 ) 

Substituting in this equation the value of/. (x) obtained from (48) we have 
^(«+l)-t{^(-)-^^n(»)}y,(») 

(49) (t = 2, .... n). 

The last line of the expression may be written in the form 

Since f. x (x), j = 1, • • •, n, is a solution of (24), this reduces to 

an expression which is identically zero. Hence in writing (49) the last line 
may be omitted. Thus we have 

(50) gi (x + l) = ±l [ f ;j (x)-j^±^ ) f 1J (x)j gj (x) (i = 2,- -.,«). 

Substituting the values of f { (x), i = 1, • • •, w, obtained from (48) in the 
first equation of system (24), we have 

(51) ? 1 (* + 1) = E^^*^*)?,(*) + E* V (*)S',(*> 
Replacing in (51) and (50) the functions yjr r (x) by their original values, we 
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obtain the transformed system of equations : 

/ l\ m i n f ( x) n 

(52) y 1 ( SB+ l)_o 1 ^l+-j g 1 (x)=Tj^ ) ^ J (x)9 l (x)+T ) ^ lj (x)g J (x), 

(53) gi {x + \)- ai (i + \Y3^) = i\h^)-j^^^^)\gM) 

m . . . ( i = 2 > ■•■•»)• 

This system has the simple solution 

i7i( a; )=/ii( a; )' ^O) = (i = 2, •••,»). 

By substitution of this solution in (48) one reproduces of course the known 
solution of (24). We shall prove next that system (52), (53) has a second 
solution, and by means of it find a second solution of (24). 

It will be noticed that the n — 1 equations (53) involve only the n — 1 func- 
tions <7,(x), i = 2, • • ., n ; and therefore (53) is a system of n — 1 equations 
in these unknowns. It satisfies the hypotheses of the theorem indicated in the 
remark at the close of § 2 ; that is, there exist a constant M' ( = M) and a D- 
region such that each of the coefficient functions 

is analytic in D and the inequality 

M'c 






■(' 



is satisfied by every x in D . The statement referring to the analytic nature of 
the coefficients is evidently true; for in a suitable Z)-region f n (x) does not 
vanish. It is clear also that there exists a constant M' such that the above 
inequality is true provided that there is a _D-region in which / ;l (x+l )//j 1 (jc+l ) 
is in absolute value equal to or less than 1. That this is true follows at once 
from inequalities (40) to (43). 

From the above considerations and the theorem indicated in the remark 
referred to may be drawn several conclusions which are important for the 
succeeding argument : 

1) There exist particular regions _P 2 , Q 2 , D 2 and constants c and M X ( = M'), 
such that inequalities analogous to (30) to (34) hold and such that also 

(54) u>2Jf 1 + l, AJv<l, |/ n (a>)// u (aO|^l, for x in i> 2 . 

2) The system of equations (53) has a particular solution, g i2 ( x ) , ?'= 2 , • • -,n y 
each function of which is analytic throughout D 2 . 

3) The increase or decrease of the functions which constitute this solution is 



118 R. D. CAKMICHAEL: difference equations [January 

bounded as follows : 

, w — 2-4, , . , , v 
alx m > -r 1 < g„(x) < a'ar"* T , 

12 ' t) — A. ^'■ /22V - ' I ] 2 '» — ^1 1 

(55) ' for a; in § 2 ; 

1^(^)1 < | fl^af | ^^ (i = 3, ...,»), 

M _ 2M — 1 , , N , , , u — 1 

l^l^*^<ly»(*>l<l^l «-jy,-i ' 

for a: in ,P„ . 

I y a («) I < I «2*™ 2 1 ^^TjfTTl (» = 3 , ■ ■ • , n) , 



(56) 



i 



These n — 1 functions are a part of a particular solution of the system (52), 
(53). To find the other function g n (x) of this solution we have the non-homo- 
geneous equation 

/ l \ »i 

( 5T ) n f (x) ^ 

= s-fS« w (*)y u (*)+ iv*)^(*)- 

We obtain first a solution of the reduced equation 

/ IX™ 1 n f (x\ 

(58) g( x +l)-a l [\+-} g^)=^ J J^ ) j> lj {x)g{x), 

an equation which satisfies the hypotheses of the theorem indicated in the 
remark at the close of § 1, as is shown by relations (54) and those indicated in 
the lines preceding (54). Accordingly it has a particular solution g(x) 
analytic and not zero throughout D 2 and approximately characterized by the 
inequalities 

I «i a*" 1 1 V ZTA < 1 9 ^ I < I a ' I ^A, ' x in & ; 

( 59 ) 

Substituting in (57) 

(60) g u (x) = h(x)g(x) 

and making use of the fact that g(x) satisfies (58), one obtains by easy reduction 

(61) A(, + 1) - h(x) = ± Y(" ) + l)^ = V(X) - 

It is clear that ri{x) is analytic throughout D 2 . 

Bounds of increase or decrease of r){x) may be readily obtained. From (55), 
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(56), and (59) we have 

i*(-)i<|(S|)' 



^7«2— 7»i —2 



(62) 



Mv . 

V-2AS fora,ln ^5 



t m%— mi— 2 



Jf(« - 1 ) . . p 



From the first of these inequalities it is clear that the series 



(63) 



«(SB) = 1j(»-l) + lj(a!-2) + 1j(!B-8) + 



is uniformly convergent throughout any closed region in Q 2 ; and hence the 
function s(x) is analytic in Q 2 . We have also 

(64) s(x+l) — s(x) = r)(x), 

so that h(x) = s(x) is a solution of (61) which is analytic in Q 2 . But its 
behavior near the positive real axis is complicated ; and the solution of (52), 
(53), and hence of (24), to which it leads is composed of functions which are 
complicated in their character in the same region. We seek therefore another 
solution, one which will avoid this difficulty near the real axis ; it will of course 
differ from s ( x ) by a periodic function. 

A solution of (61) in the form of a contour integral, the path of integration 
going to infinity along two parallel lines, will be found to serve our purpose. 
We form the integral 



where the path of integration L is yet to be chosen. Since r)(x) is analytic 
throughout D 2 , it is evident that all the infinities which the integrand has in D 2 
are at the points x ± r, where r is zero or a positive integer. Let AB be the 
straight line containing all these points. Let the path of integration L or 
CKRH lie entirely in D 2 and be formed in the following manner : It consists 



H 



R 




K 



Fig. : 



of three parts : a part RH lying entirely above AB and extending to infinity in 
a negative direction parallel to the axis of reals ; a part KG lying below AB and 
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extending to infinity in a negative direction parallel to the axis of reals ; a finite 
part KR which crosses AB once between x — 1 and x and at no other point. 
Moreover the path CKRH is such that E 2 lies entirely to the left of it. From 
(62) it is evident that the integral along such a path has a meaning. It is clear 
that for a particular x the path of integration may be deformed at will without 
changing the value of I(x), provided that during the deformation it has always 
the properties stated in the definition of CKRH. Furthermore, _Z~( x ) is analytic 
inZ> 2 . 

We show that h(x) == I(x) is a solution of (61). In view of the fact that 
the integrand considered as a function of x is periodic of period 1 it is clear that 
we have 

<w) A« + i>-Xr^£*. 

where L' is the path CTSH. Now 






But by Cauchy's theorem the last integral is equal to ??(x). Hence, subtract- 
ing (65) from (66), member by member, we have 

I(x+1)-I(x) = v (x); 
that is, 

h(x) = I(x) 

is a solution of (61). It is this particular solution which we shall take for the 
value of h(x) in the remaining discussion. It is analytic in D 2 , as we have 
seen. 

It is necessary, for later use, to know a bound to the increase or decrease of 
| h(x) | for x in Q 2 and P 2 . We consider the two regions separately, first taking 
the region Q 2 . On account of (64) and the analytic nature of both h(x) and 
s(x) in Q 2 , it is evident that h(x) — s(x) is a periodic function of x of period 
1 and analytic throughout Q 2 . Now 

(67) |A(as)|^|A(a!)-8(aj)| + |«(a!)|; 

accordingly we find upper bounds for |A(cc) — s(#)| and |s(*)| separately. 
Consider the function h(x) — s(x) for x = u± v V — 1 , x being in Q 2 , where 
v is fixed and where u runs over any stretch of length 1 . For such variation 
of a;, | A (a;) — s ( x ) | has evidently a maximum which is entirely independent of 
the range for u but depends directly on v . Call this maximum m ( ± v ) . 
Then we have 

(68) | h(x) — s(x) | = m(± v), for x in Q 2 , 

where m(± v) is a single- valued function of ± v. 
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l s ( x )l = Zlh( a; — 01 < 



Mv 



v - 2A l £- 



(sr<-" 



by (62). Hence we have 

s(^) < o-r 

1 v y| t) — 2J.j 

Therefore we obtain the result 

B, 



(a)V~|s|(5)'(.-|) 



m 2 — wij— 2 



j»2 — »»i— 2 



(69) 



! ( a; )l<i 



(5)'' 



for ar in ^ 2 , 



where B x is a properly chosen constant. Then, finally, from (67), (68), (69) 

we have 

B, 



(70) 



B \ ( a \ x 
|A(a:)| <m(±t>)+ i— 12 ( - I »"" 



, for a; in Q 2 . 



In order to find an upper bound for | h ( x ) | when a; is in _P 2 we proceed as 
follows: Let AB be the line parallel to the axis of reals and containing the 

F 
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Fig. 3. 



point x. Construct the contours ELF, EMF, ENF, P as in the figure, 
where the stretch CD is equal to J. Let us denote by I L (x), I M (x), I N (x), 
I P (x) respectively the integral 

7) ( z ) dz 



I 



1 g277t(X— Z) 

taken around these contours. For any value of a;, either I M (x) or /^(a;) repre- 
sents an integral which is convergent, though the path of integration of one of 
them may pass through an infinity of the integrand. I L { X ) an d Ip{ x ) always 
have a meaning, if it is understood that the contour P may vary so as not to 
cross AB at a point congruent to x . 
Since A ( a; ) = I L (x), we clearly have 



\h(x)\ = \I M (x)\ + Ip(x), 
\h(x)\ = \I#(x)\ + Ir{x), 



for a; in P , 
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where T P (x) is the quantity obtained by applying to I P (x) the calculus of 
residues and replacing by its absolute value each term of the sum so obtained. 
Let x vary along AB . For each x take the smaller of the quantities li^cc)! 
and |jTy(a:)|; call this quantity t(x). From the nature of the integrand in 
I(x), it is evident that t(x) has a maximum; we denote this maximum by 
H ( ± v ) , since it depends on ± v alone. Thus we obtain 

| h(x)\ = fi(d=v) + I P {x), for x in P 2 . 
Now from (62) and the calculus of residues it is clear that when x is in _P 2 we have 

M ( u -JL v I (<h \-\ x _ „ _, 

\-lh\\aJ {X l) 



\I P {x) 
Hence we obtain 



I /*»(*) I < 



u-2M, 
B. 



(;)" 



v m s - mi I forxinP,, 



where B 2 is a properly chosen constant. Therefore we have 



(70) 



\h(x)\ <p(±v) + 



B„ 



m 



for x in P. 



Since h(x) and g(x) are known functions, both analytic throughout D 2 , 
equation (60) gives us a function g 12 (x) analytic in D 2 . From (59), (70), (70) 
we have 



\g l2 (x)\<\a*x" 



(71) 



■y-m(± v) 
v — A,~ 



+ \aix m 



B., 



i / %, , ,f JL (±'v)(w— 1) 

| sr la (aj) |< | aja^ | \,_ ^_ 1 ; + | o*ai 



5. 



, for a; in Q 2 ; 



«— 1 



for x in _P 2 . 

From (48) we now have a second solution of (24) ; calling itf (x) , t = 1 , 
we may write 

/l 2 ( 3; ) = 9 , 12( a3 )' 



n, 



(72) 



f*( x ) = f$)9»( x ) + 9«( x ) 



(i = 2, •••, »). 



Also from (55), (56), (71) we have the following inequalities : 



(73) 



!A(*)I<I«^ 



5„ 



a; |' v 
v 



A, ' ' 



>(±«) 



W-^ ' 



|/ 22 (a;)|<|a^^| ?7 — ^+|/ 12 (.t)|, 



|/ M (s)| >l«;g- r t ,_ 2 ^' - l/»(s)|. 

|/ a (*)| <|aj[ a -4|--i ;r+ |/ ll («)| (j = 3, ...,n), 



for a; in Q 2 ; 



1911] AND THEIE ANALYTIC SOLUTIONS 123 

, , x , B, U—l , a(±v)(tt— 1) 



(74) 



l/„(<0 I < I «^| ,7^^-zi + \f»( x ) I ' 

for a; in P, . 



u-2M,-l 



\f 22 (x)\>\a* 2 x™> l u _ M 1 _ 1 - \f n (x)\, 



M 

\f i2 {x)\ < l^af | tf _ j^_ ! + |/»(«)1 (i=S, .-.,«), 

It should be observed that the results of this paragraph, inasmuch as they refer 
only to the region D 2 , are again independent of the nature of the functions 
^.(x)mE 2 . 

If n > 2 we may obtain a third solution of (24) in the following manner : 
The solution of (24) just obtained depends intimately upon the two smallest a's 
in the system, as inequalities (73), (74) show. Hence it follows that (53) has a 
corresponding solution depending in like manner upon the two smallest a's in 
that system ; this solution we shall call g a ( x ) , i = 2 , • • • , n . It is analytic 
throughout some region Z> 3 and satisfies the inequalities corresponding to (73) : 

, > , , , B, v , , v ■ m, ( ± v ) 

q,Jx)\ < a x »x m A 7—f, -r + a x ,x m A l -^-r- ' 

\ira\ j\^-\ 3 i up v — A„ ' 2 ' v — A. 



\g 3i {x)\<\aix m *\ v _ A +\g 23 (x)\, 



2 

V 
A 

2 



(75) for a; in Q 3 ; 

, m , ,v — ^A. , . ., 

\gn( x )\>\ a l xmz \ V _ A -\9a(«)U 



2 

A 



g i3 (x)\<\aix™°\— 2 j~+\ g2i (x)\ (i = 4, ..-,»), 



where Z) 3 , Q 3 , ^(db -u), A 2 , B 3 are defined in a manner analogous respectively 
to the definitions of D 2 , § 2 , m ( ± v ) , J. x , B l . Also there exists of course a 
set of inequalities similar to (74) for x in P 3 ; it is not necessary to write them 
out. 

We have thus n — 1 functions of a third solution of (52), (53). To find the 
other function, which we will call g l3 (x), we have an equation corresponding 
to (57) above. If we write 

g l3 (x) = h l (x)g(x), 

where g(x) is defined as before, it is easy to see that h 1 (x) is a solution of 
(76) h l (x + l)-h 1 (x) = ±H^^^ Vl (x). 
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From inequalities (59), (75) and those for x in P 3 corresponding to (75) and the 
fact that Vi( x ) * s analytic throughout D z it follows that the method used for 
solving (61) is applicable to (76) and that a solution may be obtained as before 
in the form of a contour integral. We thus find a third particular solution of 
(52), (53). Just as before this leads to a solution of (24). It is evident that 
this process may be continued until n particular solutions of (24) are obtained 
and that there exists a Z>-region in which each of the n? functions of these solu- 
tions is analytic. We shall call these n solutions 

( 77 ) flj( X )^f2j( X )^-^fnj( X ) (; = 1, •••,»). 

It is clear that we may employ (27) to extend each of these solutions across 
the plane to the left and in this way prove that each function f r (x) is analytic 
throughout the finite plane except at the singularities of ^..(x) and the points 
congruent to them on the left. 

A reference to the discussion (in the preceding section) of the bounds of 
!_£.,(!*;) | in the region E x will bring out the fact that an argument similar to 
that employed there is applicable to any one of the n particular solutions which 
we have obtained, and that the result in each case takes the same form as that 
obtained for the first solution. 

Let us now write 

f ( .(x) = a'.x^U^x) ( i,j = i, ...,„) , 

where the functions U^^x) are defined by these equations. Inequalities (40) 
to (43), (73), (74) and the analogous inequalities for the other solutions enable 
us to conclude (without further argument) that when x approaches infinity 
in the positive direction along any line parallel to the real axis we have * 

[0 if i+j, 

(78) lim £/,(*)= . . 

11 if i=j. 

This fact will enable us to prove the independence of the solutions (77). The 
necessary and sufficient condition for such independence j- is that the determinant 
of which f r (x) is the general element is not identically zero. Now we have 

\f t j ( x ) | = a\ a\ ■ ■ ■ a' n a-i+»-+-+— | U { . ( x ) \ . 

From this result and (78) it follows that the determinant of which f. .( x ) is the 

*This result brings into prominent notice the important question of the asymptotic character 
of the solutions for every approach to infinity. The full discussion of this subject requires a 
separate investigation. There are several excellent papers by Horn, Ford, Perron and others 
treating the highly important special case in which x runs over a set of points x , x +l, x -|-2, • • •. 
See references on pp. 160, 190, 191 of Horn's paper in Crelle's Journal, vol. 138 (1910). 
To the papers which he lists should be added one by Perron in Jahresbericht der deut- 
schen Mathematiker-Vereinigung, vol. 19 (1910), pp. 129-137. See also the paper by 
Watson referred to in the introduction to the present paper. 

t Cf. Pincherlb and Amaldi, Le operazione distributive, p. 218. 
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general element is not identically zero. Hence the solutions (77) are independent 
and form therefore a fundamental system of solutions of the system (24) ; and 
hence the general solution is 

n 

/<(*) = T, «>j(x)fij(x) (< = l, ■■■,»). 

where the functions a>, (x), j = 1 , • • •, n, are arbitrary periodic functions of 
period 1. 

In the theorem below, the inequalities bounding \f (x)\, i, j = 1, ■ • •, n, 
are not given ; otherwise the principal results of the section are as follows : 

The system of equations 

n n 

/«(*+l)=2>y(«).//0'0 or f i (x-l)=j:^ i .(x-l)f J (x) « = !,• ..,»), 

y{r r (x) being restricted as in the theorem of § 2, has the fundamental system 
of solutions 

fu( x )'f v ( x )' ■•■-f nj { x ) (i = l. •■■,»), 

determined in the preceding pages and having the following properties : Each 
function f,(x) is analytic throughout the finite plane except at the singulari- 
ties qfijrj.^x) and the points congruent to them on the lejt. The general solu- 
tion of the system of equations is 

n 

/«(*) = 2>,0e)/„(*) « = i,...,n), 

where the functions a> . ( x ) , j = 1 , • • • , n, are arbitrary periodic functions of 

period 1 . Further, for x approaching infinity along any line parallel to the 

positive real axis we have 

f0 if iJpj, 
lim f. . ( x ) a . x x '"■> = I 

J ' A ' J jl if i=. 



--J- 
§4. Reduction to normal form. General existence theorems. 



Let 



(79) G i (x + l) = x^a i .(x)G j (x) (i = l, •••,»), 

be a system of n first order linear homogeneous difference equations involving 
n unknown functions G 1 (x), ■ ■ •, G n ( x) of the complex variable x, the known 
quantities entering into the equation being defined as follows : 

1) a is a constant ; 

2) the functions a. . ( x ) are single- valued and 

a. J (x)'=a. j + a' ij x- l + a'ljX- 2 + ••• (»,;'= 1, •••, w) \x\> R ; 
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3) the constants a { . are such that the roots of the equation 

[OiH+j, 

(80) \ aij -S ijP \ = 0, «„=.... . 

are in absolute value different from each other and from zero. 

Equation (80) is the characteristic equation associated with (79). For con- 
venience we assume that its roots a x , • • • , a n are arranged so that 

(81) 0<|.a 1 |<|o 1 |<...<|aj. 

The object of this section is to extend the theorem of § 3 to the system in its 
more general form and then to apply the result to a single equation of the rath 
order. The work is carried out by first reducing (79) to the normal form 
already studied; In this discussion two properties of the T-function will be 
assumed ; namely, that it is analytic to the right of the axis of imaginaries, and 
that it satisfies the functional equation 

(82) T(x + \) = xT(x). 

It is obvious that a can be reduced to zero by means of the transformation 

(83) G l {x) = {Y{x)Y& i {x). 

Effecting this substitution in (79) and reducing by means of (82) we have 

(84) G i (x + l) = ±a i .(x)G } .(x) (i = i, ■••,»). 

i=i 

The first step in the normalization of (84) will be made by means of a linear 
transformation with constant coefficients. This done, a further reduction (of a 
different type) will enable us to remove the term in x~ 1 from the resulting 
coefficient functions a..(x), i 4= j. 

Let us put 

n 

(85) G i {x) = '£ t a iJ F j (x) (i = i, ••-,»), 

where the determinant | a. . | is different from zero. Making this substitution in 
(84) and solving for F i ( x + 1 ) one may write the result in the form 

(86) F i {x + l) = ita iJ (x)F.(x) (i=i, ••-,«), 

j- 1 

where the functions a i .(x'), being linear combinations of the functions a ; .(cc) 
with constant coefficients, may be expanded in powers of x~ l for | x | > It'. 

In connection with this transformation, let us consider the effect in the linear 
substitution 

n 

(87) y' i = Ha iJ y j (<=i, ••-,»), 
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of replacing y\ and y i by their values as given in the transformation 

n n 

(88) y i = H a -i J y j , y'i = T, a ijy'j (<=i, ■■■,»).. 

Having made these substitutions, solve for y\ and write the result 

n 

y"i = H%yj (i=l, •••,»). 

Comparing the operations of this paragraph with those of the preceding, it is 
easy to see that the constant term in the function a r (x) must be identical with 
a v for each pair of values i and j. But, on account of (81), it follows from the 
well-known theory* of linear substitutions that a proper determination of a., in 

(88) will transform (87) to the simpler system 

(89) £ = «!& (i = l, ■••,»). 

Suitable values of a. . may be determined * by solving the following n systems 
of equations : 

n 

(90) ll a ij * j T — a T a iT= Q (*I T = 1» •••.»). 

j=l 

a separate system being formed for each t. That each system (90) is consistent 
follows from the fact that p = a T is a solution of the characteristic equation (80). 
Hence the ratios of a. T for fixed t and i = 1 , • • • , n are determined by (90). 
Out of the infinity of proportional solutions of each of these n systems of equa- 
tions we suppose a particular one chosen, and that the quantities a., have these 
fixed values throughout the remaining discussion. 

If the values of a { . so determined are used in (85), that transformation throws 
(84) into the simpler system 

n 

(91) F i {x + \) = a.F i {x) + Y,h( x ) F i( x ) (*=l,--,»), 
where 

b iJ (x) = b. J x- 1 + b'ijX- 2 + ■■ ■ (i, j= 1, ••-, n) \x\> B. 

In order to remove the term in a; -1 from each of the coefficient functions- 
b. . ( x ) , i =|= j , we employ a transformation of the form 

(92) F i (x)=f i (x) + j:^f J (x), B u = 0, «=!,• .■,»), 

j=l x 

where the constants B { ., i =^j, are as yet undetermined. The determinant of 
this transformation is 



A = 



o« + — ' 



= 1 -f terms in x~ 2 , a; 



2 ^.-3 



* Cf. Hoen, GewShnliche Differfntialgleichungen, pp. 72-74. 
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Substituting these values for F { (x) in (91) and writing only those terms which 
do not involve x~ 2 , x~% • • • , we have, when | x | > R\ 

/,««+i)+|;(^' 1 )/ J («+i)-<.,^(«)+«,E^/ / (x) 

the accent in 2' denoting that in the summation the term for which j == i is to 
be omitted. 

It is necessary to solve (93) lor f ( (x + 1). We proceed thus : The cof actor 
A. f of the element in the ith row and jth column of A has clearly no constant 
term, except when i equals j, when there is a constant term 1 : the term in 1/x 
for i =%=j is clearly — /3 j{ /x. For, if we cross out the ith row and jth column 
of A , there remain n — 2 elements with a constant term 1 ; and therefore the 
expansion of this minor can contain only one term in l/x, namely that term 
which contains all the n — 2 factors 1 and whose remaining factor is not in the 
same rows or columns as these. Therefore this term is ± fi.Jx ; and, since an 
even number of interchanges of rows and columns bring the elements jS../x and 
Pjjx to the place in the first row second column and second row first column 
respectively, the negative sign must be chosen. In the cofactor A tt there is no 
term in l/x; for every term in the expansion containing one factor l/x neces- 
sarily contains another. Hence the cofactors may be written 

/3.. 
A.j = s - -f terms in x~ 2 , ar 3 , • • • , ( i + j ) ; 

A w = 1 + terms in x~ 2 , a; -3 , • • • . 
Solving (93) for f.(x + 1 ) we readily obtain 

(94)fr,(»+l)-(«, + ^ + ---)^^ 



A denoting the determinant 



(» = !,■■-, n), 



A= S .+ 



x + 1 



For values of x such that | x \ = R" , It" being a properly chosen constant greater 
than R' and greater than 1, the determinant A is everywhere different from zero 
-and is expansible in the form 1 -f- terms in x~ 2 , x~ 3 , • • • . If we divide by 
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A we obtain therefore 

| x | = It : 

the terms of order x~ 2 , a; -3 , • • • in the coefficients of /.(as) in the right hand 
member are of course different from the corresponding terms in (94). 
On account of (81) it is always possible to determine /3 r , i 4=i, so that 

(fl 4 -o,)^ + 6 # = 0. 

When the values of /3 r so determined are substituted in (95), that system reduces 
to the form 

(96) ^+ i )=(««+^)/ 1 (»')+i:(^+---)/,(«) (i=i, ...,») 

|a?|S2?". 
Now 



/ 1\*« &.. 

z. ( 1 + - )"■ = a. + — + 
' \ a;/ * a; 



and hence, changing to more convenient notation and employing relations (81), 
we obtain the normal system 

(97) / l ( a> + l)_ a( ^l + _j / 4 (a:)_2:^(*)/,(*) (*=l,-,-), 

where 

o<l«il <!«■! < ••■ <KI> m i = ^> 

and 

tf>..(as) = 0;.x- 2 + ^>- 3 +... |*|==i?, 

where H = 12". This is the normal form which we have studied in the preced- 
ing sections. 

Combining the transformations (83), (85), (92) we readily obtain 



(98) G ( (x)={T(x)}-± « a [/>)+£^ /,(*)] 



(i = l, •••, n). 



It is evident therefore that if we replace / { ( as), i = 1, •••,«, by the particular 
solution/^. (as), i = 1, • • ., w,^' fixed, we obtain a particular solution of (79) 
which we will call G v (as), i = 1, • • •, «. Doing this for J = 1 , • • ., n, we 
have w particular solutions 

Gy(x), G 2J (x), -.., £„(») U=i,-.»). 

Trans. Am. Math. Soc. 9 
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where 

m r n 8 ~i 

(99) G ij (x)={T(x)}«Z« ik [f k .(x) + j:^f lj (x)\ «,i = i, •••,»). 
Let us associate with (79) the corresponding equation 

(100) G i (x-l) = ±A ij (x-l)G J (x) (.■=!,. ..,.). 

Now there is a .P-region, as we have seen, in which the functions fij(x) are all 
analytic. It follows from (99) that the functions G tJ (x) likewise are all 
analytic in the same region. By means of (100) the functions of each of these 
solutions may be extended across the plane to the left ; and in this way it may 
be shown that each of the n 2 functions G (j (x) is analytic throughout the entire 
plane except at the singularities of A ij (x) and the points congruent to them on 
the left. 

Let us again write 

/«(*) = a^'U^x) (i,i = l, •■-.»). 

If we substitute these values for/^(») in (99) we have 

G ij {x)={T{x)Ya^V tj {x), «,j=h- ■■,«), 



where 



^(*)=g««[^(*) + i:§- , ^(*)]. 



From (78) it follows that when x approaches infinity in the positive direction 
along any line parallel to the real axis we have 

WmV ij {x) = a ij . 

This result and the fact that the determinant | a v | is different from zero lead 
easily to the conclusion that the n particular solutions G tj (x) are independent. 
The method of proof is that employed for showing the independence of the solu- 
tions/^ (x) in the preceding section. Hence the general solution of (79) is 

n 

G t (x) = Y.°>j(x)G ! ij(x) (i=i, ••■,»), 

where the function a>j ( x ) , j = 1 , • . • , n , are arbitrary periodic functions of 
period 1. 

The principal results of the section may be stated in the following theorem : 

Let there be given the system of equations 



G i (x+l) = x"J2a ij (x)G J (x) (» = i, 

1) a is a constant ; 



where 
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2) the functions a tJ (x) are single- valued and 

a {J (x) = a t . + a'^.x- 1 + a'!.x-' 1 + ■ ■ ■ (i,j= 1, • • ■, ra) |aj| > R'; 

3) the constants a t . are such that the roots p = a. of the equation 

\a.. — S..p| = 

are in absolute value different from each other and from zero. 
Associate with the given system the corresponding one 

n 

G i (x-l) = *£A iJ (x-l)G J (x) « = i, ■-,«). 

Then the given set of equations has the fundamental system of solutions 

#«(*)> G 2J (x), -.., G nJ (x) (i = l, ■■-,«), 

determined in the preceding pages and having the following properties : Each 
of the n 2 functions is analytic throughout the finite plane except at the singu- 
larities of A..(x) and the points congruent to them on the left. The general 
solution of the system is 

n 

G i {x)=^ t co.(x)G iJ (x) (i = i, ••-,»), 

where the functions w.(x),j = 1, • ■ ■, n, are arbitrary periodic functions of 
period 1. Further, if x approaches infinity in the positive direction along any 
line parallel to the real axis we have 

Km G.. ( x ) { r ( x ) } - aj' x~ m ' = a.. , 

where the constants a tf are determined by equations (90) and the constants m } 
are those introduced in (97). 

Let us now consider the single difference equation of the nth order 

H(x + n) + (x + n) a b l {x)H{x + n— 1) + (a? -f n) a (x + n — l) a 
(101) 

xb 2 (x)H~(x + n — 2) -\ f- (x + n) a - ■ ■(x + l) a x a b n (x)JI(x)=0 

where a is a constant and b f (x) is single- valued and 

b.(x) = b i +b' i x- x + b';x- 2 + ••• (» = 1, ..., n) |*| > 5, 

and where the constants b { are such that the roots a l , ■ ■ •, a n of the equation 

(102) p» + 6,,- 1 + 6 2/J »- 2 + • • • + b n _ lP + 6 n = 

are in absolute value different from each other and from zero. 
Making the substitution 

(103) H(x) = {T(x) yfi(x) 
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we have 

(104) H(x + n) + \(x)H(x + n-\) + ... + b n (x)H{x) = 0. 

It is evident that this equation is equivalent to the following system 

n 

^{x+l) = ^{-b j {x)}H J {x) = ^ 

(105) j=i 

H i (x + \)=H i _ l (x) (i = 2, ••,»), 

where 

(106) H n (x) = H(x). 

Now the characteristic equation [compare (80)] associated with (105) is readily 
reduced to the form (102). 

If we solve (105) for H t (x) in terms of^(cc + l),i=l, • • -,n, and in the 
result write x — 1 for x , we have 

H i {x-l)=H i+l (x) (i = l,...,B-l), 

Consequently, from the foregoing theorem, we have the following conclusions : 
The system (105) has a fundamental system of solutions 

H v {x), H 2J (x), ...,H n .{x) o = i,. ..,«), 

each function of which is analytic throughout the finite plane except at the sin- 
gularities of the functions 

KW KW ° ' "' ,n) ' 

and points congruent to them on the left ; further, if x approaches infinity in the 
positive direction along any line parallel to the real axis, we have 

lim If Ax) aj x x~ m > = a. . , 

where a tj and m i are constants corresponding to the a tJ and m j of the above 
theorem. Since only the ratios of cL , j fixed, i = 1 , • • • , n , are determined, 
we take a nj , = 1 , j = 1 , • • • , n . 

Combining (103) and (106) we have 

H{x)={Y(x)}*H n {x) 

Hence we have a fundamental system of solutions Il.(x),j = 1 , ■ • •, n, of (101) 
in the form 

mx)={r(x)}°ii nj (x) o-=i, ...,„). 
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Making use of this equation and the properties of H nJ (x) stated above and 
determining the position of the singularities from the equation itself by extend- 
ing the solutions across the plane to the left, we obtain the results stated in the 
conclusion of the following theorem : 
Given, the difference equation 

H(x+n) + {x+nyb l (x)H(x + n-\)+ ■ ■ ■ + {x + ny- • .(x+iyx°b n (x)H(x)=0 

or 

PJ,x)H(x+r,) + P n _ l {x)H{x+n-\)+..-+Plx)H(x+\)+H(x) = <S, 

where a is a constant and b t (x) is single-valued and 

0i (x) = b t + fyx- 1 + %x~ 2 + ... (»=l,...,n) |os|>i?, 

and where the constants b t are such that the roots a 1 , ■ ■ •, a n of the equation 

p' + & lA -i+... +b n _ lP + Z>„ = 

are in absolute value different from each other and from zero. 
This equation has the fundamental system of solutions 

H,{x), H 2 {x), ..-, H n {x) 

determined above and possessing the following properties : Each of the func- 
tions /7,(o5), • • -, II n (x) is analytic throughout the finite plane except at the 
singularities of the functions /S i (x),i= 1 , •••,«., and the points congruent to 
them on the left. If x approaches infinity in the positive direction along any 
line parallel to the real axis, we have 

lira II. (x){T(x)}- a aj x x~ m '= 1 {j = l,---,n), 

where m, are constants. 

The general solution of the eqation is 

II{x) = co^H^x) + •■• + w n {x)H n {x) 

where ^(a;), • • •, a> n (a;) are arbitrary periodic, functions of period 1. 

It may be observed that a modification of the foregoing method will lead to a 
different fundamental system of solutions of (79) and likewise of (101). We 
confine our attention to (79). A few remarks relative to the equation 

g(x + l)-a\l + -J g(x) = v (x) 
will sufficiently indicate the necessary changes. Employing the notation 
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we have two formal solutions of the above equation as follows : 

, x f- v(x + i) 

CO 

g(x) = T)(x— l)+^X(x- l)\(x — 2)-- -X(x — i)r)(x — i — 1). 

The first of these, as we have seen, leads to a solution of (1). The second may 
also be employed in a similar way to obtain a solution of the same equation. 
The extension to a system of n equations is obvious. In this way a second 
fundamental system of solutions of (79) may be found. 

Suppose further that transformation (83) is modified by the use of a func- 
tion r ( x ) which satisfies the equation 

T(x + l) = xT(x) 

and is analytic to the left of the axis of imaginaries, this function replacing 

T(a;) in (83). This modification having been made, it will turn out that each 

function of the fundamental system just found is analytic throughout the 

finite plane except at the singularities of a..(x) and the points congruent to 

them on the right ; and their properties may be developed by a method similar 

throughout to that employed in the paper and with like results. 

Princeton University, 
August, 1910. 



